Simulation study of BESIII with stitched CMOS pixel detector using ACTS* 
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Reconstruction of tracks of charged particles with high precision is very crucial for HEP experiments to 
achieve their physics goals. As the tracking detector of BESIII experiment, the BESIII drift chamber has suffered 
from aging effects resulting in degraded tracking performance after operation for about 15 years. To preserve 
and enhance the tracking performance of BESIII, one of the proposals is to add one layer of thin CMOS pixel 
sensor in cylindrical shape based on the state-of-the-art stitching technology, between the beam pipe and the drift 
chamber. The improvement of tracking performance of BESIII with such an additional pixel detector compared 
to that with only the existing drift chamber is studied using the modern common tracking software ACTS, which 
provides a set of detector-agnostic and highly performant tracking algorithms that have demonstrated promising 
performance for a few high energy physics and nuclear physics experiments. 
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I. INTRODUCTION 


The Beijing Spectrometer (BESIID) [1] at the Beijing 
Electron-Positron (BEPCII) Collider has been a great success 
producing prosperous physics results [2] in the 7-charm sec- 
tor since 2009. However, after operation for about 15 years, 
the detectors at BESIII have been subjected to aging effects, 
which result in degraded performance of the detectors [3]. 

The tracking detector at BESIII is the Multilayer Drift 
Chamber (MDC), which provides measurement of momen- 
tum and position of the charged tracks, and information of 
energy loss in unit path length, i.e. dE/dx [4], of the charged 
tracks for particle identification. As shown in Ref. [3, 5], 
due to the beam-induced background with a hit rate up to 
2 kHz/cm?, the gain of the MDC cells in the first ten lay- 
ers has shown an obvious decrease with a maximum decrease 
of about 39% for the innermost layer cells in 2017, which 
further leads to reduction of the spatial and momentum res- 
olution of the charged tracks. Since BESIII is not expected 
to complete its mission in foreseen years [2, 3], the track re- 
construction performance must be preserved and enhanced in 
order to not compromise the physics goals of BESIII, along 
with possible upgrade of the tracking system with state-of- 
the-art technologies for detection of charged particles. To get 
well prepared for the potential malfunction of the MDC due 
to aforementioned aging problem, plans of upgrading the BE- 
SIII inner tracker based on different technologies have been 
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proposed [3]. This includes replacement of the inner tracker 
with a new inner drift chamber [6] or a cylindrical gas elec- 
tron multiplier (CGEM) tracker [7], which has attractive fea- 
tures such as high counting rate tolerance and low sensitivity 
to aging. 


Compared to the gaseous detectors, the silicon pixel detec- 
tor has excellent spatial resolution and good radiation resis- 
tance. Therefore, one of the additional options for the BESIII 
inner tracker upgrade is to use large-area thin Complementary 
Metal Oxide Semiconductor (CMOS) pixel sensor with good 
spatial resolution based on the cutting-edge stitching technol- 
ogy, which has already been used to produce CMOS pixel 
sensors for medical imaging applications [8-10]. Recently, 
the studies of designing a first wafer-scale stitched sensor pro- 
totype, the MOSS (Monolithic Stitched Sensor) chip, towards 
an improved vertex detector, i.e. the ITS3, at ALICE experi- 
ment, are presented in Ref. [11]. 


A Common Tracking Software (ACTS) [12] is à common 
High Energy Physics (HEP) software, providing a set of 
detector-agnostic, performant and modular tools for the track 
reconstruction in HEP using modern software technologies to 
facilitate concurrency, usability, maintenance and extendabil- 
ity to tackle the tracking challenges foreseen in HEP in future. 
So far, ACTS has been used for the track reconstruction at AT- 
LAS [13], FASER [14], sPHENIX [15], STCF [16] etc., for 
different types of tracking detectors. In particular, promising 
tracking performance of ACTS for a tracking system with a 
drift chamber is firstly presented in Ref. [16]. 


In this study, the tracking performance of the BESIII MDC 
with an additional one-layer stitched cylindrical CMOS pixel 
detector inserted between the beam pipe and the inner wall of 
the MDC is studied using ACTS as the tracking software. The 


performance is compared to that of the current BESIII track- 
ing detector with only the MDC based on the BESIII Offline 
Software System (BOSS) [17]. The manuscript is organized 
as follows. In Section II, the BESIII MDC and the proposed 
pixel detector based on the cylindrical CMOS pixel sensor are 
introduced. Section III introduces the tracking strategies in 
BOSS and ACTS. Improvement of the tracking performance at 
BESIII with an additional pixel detector is presented in Sec- 
tion IV. A brief conclusion is provided in Section V. 


Il. THE BESIII TRACKING DETECTOR AND STITCHED 
CMOS PIXEL DETECTOR 


The BESIII MDC is a cylindrical chamber operating with a 
helium-based mixture gas (He/C3Hg = 60:40) and immersed 
in a 1 T magnetic field. The inner radius and outer radius of 
the MDC is 64 mm and 819 mm, respectively. The length 
of the wires ranges from 774 mm for the innermost layer to 
2400 mm for the outermost layer. The drift cells which are 
almost square shaped are arranged in 43 circular layers, which 
alternate between the stereo layers and axial layers, i.e. in the 
order of 8 stereo layers, 12 axial layers, 16 stereo layers and 7 
axial layers. The cell dimension is about 12 mm x 12 mm for 
the 8 inner layers and 16.2 mm x 16.2 mm for the 35 outer 
layers. The inner chamber is composed of the 8 inner layers 
and the outer chamber is composed of the 35 outer layers. 

The CMOS pixel sensor has been widely used for vertex 
detectors in HEP, due to its excellent spatial resolution down 
to a few um, tolerance of high hit rate up to 108 Hz/cm?, 
and good detection efficiency and radiation resistance. For 
example, the CMOS pixel sensor has been used for the inner 
tracker or vertex detector of the STAR experiment [18], the 
ALICE experiment [19] and the sPHENIX experiment [20]. 
It's also considered as a primary choice for the pixel sensor 
technology for the CEPC physics program [21]. Despite all 
the advantages of the CMOS pixel sensor, the tracking resolu- 
tion of a traditional CMOS pixel detector, in particular for low 
momentum tracks, can be limited by the irreducible material 
budget of the support structure and cooling pipe, as shown 
in Ref. [22], where the MAPS-based CMOS pixel sensor is 
proposed for the upgrade of BESIII inner chamber. There- 
fore, tracking based on traditional CMOS pixel sensors can 
be challenging at BESIIL, where tracks with transverse mo- 
mentum down to 150 MeV are required to be reconstructed 
with good resolution. 

In recent years, the innovative stitching technology 
emerged with the aim of producing large-area sensor in one 
wafer, which can be further thinned to about 50 jum. With 
such thickness, the wafer can be curved into cylindrical shape, 
hence vastly simplifying the support structure and reducing 
the material budget. Based on such stitching technology, it's 
possible to construct a vertex detector using large-area thin 
CMOS pixel sensors with cylindrical shape and simplified 
support structure. As discussed in Ref. [11], the material bud- 
get per pixel layer is foreseen to be reduced to 0.05% Xo, to 
allow bending of the silicon. In this study, a vertex detec- 
tor which consists of one layer of stitched cylindrical pixel 


sensor with thickness of 50 jum (0.05% Xo), as proposed 
in Ref. [11], and a length of 220 mm (including insensitive 
area within 7 mm to each end), resolution of 8.66 um in the 
r-@ direction and 57.74 um in the z direction, inserted be- 
tween the beam pipe, which has a diameter of 63 mm [1], and 
the MDC inner wall, is considered. The stitching technol- 
ogy foresees cooling by air flow and no readout circuits in the 
active area. Therefore, the additional material budget from 
the cooling and readout infrastructure as well as the support 
structure in the active region of the pixel layer is expected 
to have a negligible impact on the tracking performance, and 
hence not taken into account in this study. 

The geometry of the BESIII MDC and the pixel detector is 
shown in Figure 1. 


Fig. 1. A clipped view of geometry of the BESIII MDC (cell wires 
are shown in green colors) and the pixel detector (shown in orange 
color) outside of the beampipe. 


IIl. TRACKING WITH BESIII OFFLINE SOFTWARE 
AND ACTS 


The simulation study is based on samples generated using 
BOSS, which is the offline software framework for event gen- 
eration, simulation, reconstruction, performance validation, 
physics analysis and visualization [23] for the BESIII experi- 
ment. 

In BOSS, the production of both charmonium resonances 
and continuum process from the e*e~ collisions is provided 
by the KKMC [24] generator and the decays of particles are 
modelled with EVTGEN [25], which can also be used to gener- 
ate single particle samples. The detector description is based 
on the Geometry Description Markup Language (GDML) [26] 
and the interaction of particles with detectors is simulated 
with GEANTA [27]. 

The performance of the BESIII MDC with the additional 
pixel detector is studied using the ACTS software, and com- 
pared to the performance of the BESIII MDC studied using 
the BESIII tracking software within BOSS. The workflow for 


studies of the tracking performance is illustrated in Figure 2. 
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Fig. 2. The workflow of tracking performance studies based on BOSS 
and ACTS. 


A. Tracking in BOSS 


The trajectory of a charged particle in magnetic field is pa- 
rameterized using the helix track parameters at BESIII, as 
described in Ref. [28]. Track finding and track fitting are 
two tasks in BOSS track reconstruction. Track finding is 
a pattern recognition problem of classifying measurements 
into subsets and creating track candidates, while track fit- 
ting is an estimation of the helix parameters. Two basic track 
finding algorithms are implemented in BOSS: the template 
matching algorithm (PAT) [29] and track segment finder al- 
gorithm (TSF) [28]. Both a specialized track finding method 
called TCurlFinder [30] and the global track finding method 
based on Hough transform (HOUGH) [31] have been im- 
plemented to salvage low transverse momentum tracks with 
pr < 0.12 GeV. Besides, an extended segment construction 
scheme which can achieve higher efficiency for low trans- 
verse momentum tracks has been developed [32]. The found 
tracks using these algorithms are combined and fed into the 
Kalman Filter [33] algorithm for track fitting. Further track 
extrapolation is performed to obtain the track parameters at 
other subdetectors at BESIII. 


B. Tracking with ACTS 


A detailed introduction to ACTS, including geometry de- 
scription, parameterization of track parameters and measure- 
ments, and tracking algorithms etc., can be found in Ref. [12]. 
The implementation of ACTS for track reconstruction for the 
BESIII tracking detector is outlined below. 

To use ACTS for track reconstruction, the detector geom- 
etry with detailed description of detector material and place- 
ment, e.g. the GEANT4-based detector geometry, needs to be 


transformed into ACTS internal geometry, i.e. ACTS tracking 
geometry, which has a simplified description of the passive 
detector material to facilitate fast navigation and track recon- 
struction. The TGEO [34] plugin in ACTS is used to transform 
the TGEO version of the BESIII detector geometry, which is 
created based on GDML exported from the GEANT4-based ge- 
ometry, into ACTS tracking geometry with material mapped to 
the auxiliary surfaces. The one pixel layer of the pixel detec- 
tor is converted to a cylinder surface in ACTS and the 43 layers 
of signal wires of the MDC are converted to 43 layers of line 
surfaces, dedicated to representation of measurements for a 
drift chamber or track parameters at the interaction point, in 
ACTS. 

The hits on the MDC after digitization are transformed into 
one-dimensional measurements associated to the line surfaces 
in ACTS describing the drift distance of a drift chamber. The 
hits on the pixel detector after simulation are transformed into 
two-dimensional measurements by smearing the simulated 
hits with the resolution of the pixel detector using Gaussian 
functions. Those measurements are associated to the ACTS 
cylinder surfaces, where the local x (y) coordinate represents 
the r-ọ (z) in the cylinder frame, i.e. r represents the radius of 
the cylinder, y represents the azimuthal angle of the position 
on the cylinder and z is the coordinate in the z direction. 

The BESIII magnetic field is transformed into ACTS inter- 
polated magnetic field using the field map of BESIII magnetic 
field. With an interpolated field provider in ACTS, the value 
of the magnetic field for any given position is calculated by 
interpolating from a grid of known values, e.g. eight corner 
points of a field cell in three-dimensional coordinate system. 

The track fitting and finding is performed simultaneously 
using the combinatorial Kalman Filter [35] algorithm in 
ACTS, which has the capability of rejecting noise hits based 
on the x? calculated using the distance between the hit and 
the predicted track parameters, and the covariance of the hit 
and the predicted track parameters. If multiple hits are found 
to be compatible with the predicted track parameters, the hit 
with the best x? is used to filter the track parameters. The 
smoothed track parameters at the first measurement are even- 
tually extrapolated to the beam line to obtain the estimated 
track parameters at the interaction point. 


IV. TRACKING PERFORMANCE STUDIES 
A. Monte-Carlo sample generation 


The impact of the pixel detector on the resolution of the 
track parameters is studied using single jz~ and single m~ 
samples. The samples are generated with fixed transverse 
momentum pr, cos (0 is the polar angle) uniformly dis- 
tributed between [-0.8, 0.8] , and azimuthal angle à uni- 
formly distributed in the range of [0, 27]. The impact of 


3 Suppose the wire has a direction of w and the track direction is t in the 
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the pixel detector on the tracking efficiency is studied using 
(3686) — ata Jf, J/ — pt pom events generated at 
y's = 3.686 GeV. The two-dimensional distributions of cosó 
and pr for u and 7 in this process are shown in Figure 3. 

The random background hits in the MDC due to beam re- 
lated background or electronic noises are generated in the 
standard simulation sample production at BESIIL. Since ded- 
icated background noise model for the pixel detector is not 
available yet, no background noise is considered for the in- 
serted pixel layer in this study. 


B. Track reconstruction performance 


The resolution of the impact track parameters, do, d; and 
the relative resolution of the transverse momentum pr as a 
function of pr for the BESIII MDC with an additional pixel 
layer, which is placed with the radius rpixe1 at three different 
values, i.e. 35 mm, 45 mm or 55 mm, compared to those of 
BESIII MDC, are shown in Figure 4. While ry;4 has little 
impact on the resolution of pr, the spatial resolution is better 
with smaller rpixe. With rpixe) = 35 mm, the resolution of do 
and d, is 60 um and 120 um for u` and m~ with pr = 1 
GeV, respectively, and the relative resolution of 44% for pr 
at pr = 1 GeV is achieved. The resolution of do, dz, pr 
can be improved by up to 67%, 93% and 32%, respectively, 
if the proposed Pixel detector is added to the current BESIII 
tracking detector. 

The comparison of the tracking efficiency, which is defined 
as the ratio of generated charged particles which have cor- 
responding reconstructed track in the detector region (lcos6l 
« 0.93), between the BESIII MDC with the additional pixel 
layer and the BESIII MDC only is shown in Figure 5. The ad- 
ditional pixel layer can significantly improve the tracking ef- 
ficiency for tracks with in particular low momentum or large 
Icos@l. For u with pr below 0.6 GeV in the (3686) > 
vm J/w,J/y — wtp process, the improvement can 
reach 9%. The improvement can reach 10% for m with 0.85 
< Icos@l < 0.93. The pixel with ry; = 35 mm provides better 
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tracking efficiency than other two choices. 


V. CONCLUSION 


After operation for about 15 years, the tracking detector 
of BESIII experiment has suffered from aging effects and 
needs to be upgraded to preserve its tracking performance for 
BESIII to fullfill its remaining physics goals in the next few 
years. The possibility of inserting a one-layer pixel detector 
using the large-area thin CMOS pixel sensor in cylindrical 
shape based on the stitching technology, between the beam 
pipe and the MDC, was considered. The spatial and momen- 
tum resolution was studied using ACTS and compared to the 
performance of the current BESIII tracking detector obtained 
using the BESIII offline software. The proposal based on the 
stitched cylindrical CMOS pixel sensor was found to be very 
promising, i.e. the resolution of do, d;, pr can be improved 
by up to 67%, 93% and 32%, respectively, and the tracking 
efficiency for particles in the (3686) — m ^m J/v, J/v > 
u^. events can be improved by up to 10% in the low pr 
or large lcos6l region, with the additional pixel layer placed at 
an radius of 35 mm. Meanwhile, it's the first time that ACTS 
was used for track reconstruction for a drift chamber with cal- 
ibrated measurements and realistic noise hits. In future stud- 
ies, the noise model of this state-of-the-art pixel detector will 
be studied and the tracking performance taking into account 
the noise on the pixel detector will be investigated. 
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Fig. 5. The tracking efficiency of p (top), 7 (bottom) in the process of 7)(3686) — ata J/w, J/p — mutu” as a function of particle pr 
(left) and cos@ (right) for the BESIII MDC only (black circle), and BESI MDC with an additional pixel layer (denoted as "Pixel") placed 
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